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Chapter  "I 
IHTRODUCTIOIT 


Statement  of  the  Prohlem 

Air  Force  weapons  systems  recei-ve  their  material 
support  from  five  Air  Logistics  Centers  (ALC)  operated  by 
the  Air  Force  Logistics  Command  (AFLC).  These  ALCs  are 
linked  together  and  to  ma;]or  Continental  United  States 
(COHUS)  installations  through  various  contract,  common 
carrier,  and  government  transportation  systems.  The  primary 
freight  network,  LOGAIE,  is  a contract  air  cargo  system 
which  will  cost  the  Air  Force  49.8  million  dollars  to 
operate  in  fiscal  year  1977  alone  (7)« 

The  LOGAIH  route  structure  consists  of  a primary 
(trunk)  network  which  connects  the  ALCs  to  the  Aerial  Ports 
of  Embarkation  (APOEs)  and  several  secondary  (feeder)  net- 
works which  connect  the  ALCs  to  the  customer  bases 
(18;A-2-A-4).  The  entire  structure  is  presented  in  Figure 
1 . Daily  service  is  provided  to  each  of  these  customer 
bases  via  a feeder  route  which  originates  and  terminates  at 
an  ALC. 

Attempts  to  optimize  the  operation  of  this  network 
over  the  past  decade  have  resulted  in  only  partial  success 
due  to  the  size  and  complexity  of  the  overall  system.  The 
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Figure  1 

LOGA-Ifi  Route  Structure 


j 
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LOGAIR  feeder  route  structxire  is  dereloped  by  a manual 

process  even  tbough  the  primary  trunk  network  has  been 

significantly  improved  through  time-saving  computer 

assisted  models  (11  il).  Since  this  manual  pTOcess  depends 

upon  heuristics , the  feeder  route  structures  which  are 

-1 

developed  are  not  necessarily  optimal  in  terms  of  either 

effectiveness  or  efficiency  and  consequently  the  process 

2 

is  not  conducive  to  sensitivity  analysis. 

This  thesis  develops  an  optimizing  method  of  deter- 
mining LOGAIE  feeder  route  structures  that  allows  for  a 
general  comparison  of  distribution  options  such  as  mixing 
transportation  modes  and  varying  performance  levels. 

Background 

LOQAIS  is  the  U.S.  Air  Force  contracted  commercial 
air  carrier  system  which  provides  daily  delivery  of 
reparable  and  consumable  items  to  all  major  USAF  instal- 
lations in  the  COHUS.  The  system  provides  daily  scheduled 
service  among  60  stations  including  the  ALCs , customer 
bases,  and  APOEs. 

'I 

An  optimal  solution  is  a feasible  solution  that 
is  the  most  favorable  in  terms  of  cost  (8:32). 

2 

Sensitivity  analysis  means  to  vary  one  or  more 
parameters  over  some  inter'/al(s)  to  see  when  the  optimal 
solution  changes  (8:195). 
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Air  Force  Manual  76-1 , The  LOGAIR  Traffic  Manual, 
lists  the  ob;]ectives  of  the  system  as  follows; 

1.  Establish  and  maintain  a cargo  airlift  service, 

2.  Improve  the  timeliness  and  effectiveness  of 
logistical  support  by  expanding  and  improving  the  utiliza- 
tion of  air  transportation,  and 

3.  Improve  the  quality  and  reliability  of  the 
system  (18:3-1). 

The  LOGAIE  system  supports  the  movement  of  all 
priority  1 and  priority  2 cargo. ^ Priority  3 cargo  may  also 
be  moved  to  the  first  downline  station  on  a space  available 
basis  (18:^1).  The  system  is  completely  airlift  dedicated 
and  does  not  consider  alternate  modes  of  transportation. 

Developing  requirements  for  airlift.  Each  major  Air  Force 
installation  submits  a forecast  of  priority  1 and  priority 
2 tonaage  based  upon  first  quarter  historical  data  supplied 
by  AFLC.  The  LOGAIE  Requirements  Branch  then  develops, 
from  these  data,  a series  of  requirements  matrices  between 


•^Transportation  priori-fcy  codes  are  determined  from 
a matrix  of  Force/Activity  Designator  Codes  versus  Urgency 
of  Need  Codes.  Resulting  combinations  are  prioritized  as 
transportation  priorities  1,  2,  and  3.  (Source;  DoD 
Directive  4410.6) 

4 

LOGAIE  Requirements  Branch,  AFLC/liOTSL, 

Wright -Patters on  AIB  OH. 
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origin  stations  and  destination  stations.  The  matrices  are 
then  partitioned  into  trunk  stations  and  feeder  stations."^ 

A route  structure  is  then  developed  which  considers  the 
tonnage  and  type  of  aircraft  necessary,  delivery  times, 
number  of  landings , and  flight  scheduling.  The  route  struc- 
ture is  costed  in  accordance  with  Civil  Aeronautics  Board 
(CAB)  rates  and  anticipated  contract  costs,  and  then  sent 
to  the  Air  Staff  for  review  and  approval  (13:0-2). 

Other  attempts.  Research  efforts  in  this  area  have  provided 
both  relief  and  insight  into  the  LOGAIE  problem.  Captain 
Michael  P.  McPherson  and  Captain  Brian  O'Hara,  in  a 
master's  thesis  presented  to  the  Air  Porce  Institute  of 
Technology  in  June  1976,  expanded  upon  a previous  master's 
thesis  (14)  in  developing  a computer-assisted  method  for 
deteimining  LOGAIE  route  structures.  These  studies  only 
considered  the  trunk  route  structures  in  developing  a 
mixed- integer  programming  model  for  minimizing  costs.  The 
McPherson  and  O'Hara  model  showed  that  successful  reduction 
of  LOGAIE  contract  costs  was  possible.  However,  even  though 
only  12  trunk  nodes  were  considered,  the  computer  algorithm 

5 

''^In  general,  a trunk  station  is  one  which  processes 
approximately  2,000  tons  or  more  of  cargo  per  year.  A 
feeder  station  processes  less  than  2,000  tons  (7). 
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into  two  sub-problems  (11:25-27).  This  limitation  high- 
lights the  enormity  of  the  task  of  developing  an  overall 
optimizing  solution.  Due  to  present  state-of-the-art  and 
computer  resource  constraints , obtaining  a total  optimum 
solution  for  all  feeder  nodes  in  consonance  with  the  trunk 
structure  appears  impossible. 

A study  conducted  in  1976  for  the  Department  of 
Defense  by  E.  A.  Narragon  and  J.  M.  Neil , Logistics 
Management  Institute,  contained  two  significant  innovations: 
the  study  was  (1)  the  first  formal  attempt  to  evaluate  the 
savings  possible  by  modifying  the  present  feeder  network 
structure,  and  (2)  the  first  study  to  consider  the  time/cost 
trade-offs  possible  by  utilizing  a mixture  of  air  and  sur- 
face (truck)  transportation.  The  Logistics  Management 
Institute  study  pointed  out,  "There  is  an  excellent  surface 
transportation  system  within  CONUS  upon  which  LOGALR 
managers  are  not  capitalizing  [15:G-^]."  Considerable 
savings  might  be  realized  from  an  optimum  mix  of  air  and 
surface  transportation  which  can  maintain  the  required 
performance  level. ^ 


Performance  level  is  measured  with  respect  to  time 
and  consistency  (1:28).  The  performance  level  required  for 
priority  1 and  priority  2 cargo  is  one  day  delivery  (DoD 
Directive  4410.6). 
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Justification 


Alth-ough  lowest  total  cost  expenditure  was  not  an. 
explicit  objective  of  the  LOGAIR  charter  circa  1969 
(18:3-1),  justification  for  reducing  LOGAIR  costs'"  has  sub- 
sequently been  documented  (11:6).  Soaring  fuel  costs  and 
dwindling  fuel  reserves  dictate  that  the  system  be  con- 
stantly scrutinized  for  ways  to  increase  the  overall 
efficiency  in  terms  of  defense  dollars.  The  fiscal  year 
'1968  Department  of  Defense  budget,  which  was  being  formu- 
lated at  the  same  time  the  LOGAIR  manual  was  being  written, 
comprised  45.6  percent  of  the  Federal  Budget.  In  contrast, 
President  Ford’s  fiscal  year  1977  budget  allocated  only 
25.4  percent  to  the  military.  The  Air  Force's  share  of  the 
Department  of  Defense  budget  decreased  from  55.0  percent  to 
27.7  percent  in  the  same  period  (19:108). 

Previous  efforts  to  reduce  distribution  costs , with 
respect  to  LOGAIR,  have  been  directed  at  the  design  of  the 
trunk  route  structure  in  order  to  provide  a more  efficient 
solution  by  reducing  the  total  air  miles  flown  (11:7).  No 
progress  has  been  made  in  terms  of  viewing  the  overall 
savings  possible  if  LOGAIR.  is  viewed  as  a complete  system 
itself. 
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'^The  most  recent  cost  figures  show  that  the  total 
contract  cost,  including  fuel  subsidies,  for  fiscal  year 

1976  was  $46.4  million.  The  estimated  cost  for  fiscal  yea 

1977  (disallowing  increased  fuel  costs)  is  $49.8  million 
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The  capability  to  analyze  alteinatives  within  the 
ATTjC  distribution  system,  including  various  mixes  of  air 

and  land  transportation,  inventory  costs  versus  pipeline  ^ 

times,  transshipment  points,  etc.,  is  necessary  if  the  most 
cost  effective  system  at  any  one  time  is  to  be  identified. 

The  performance  of  the  total  distribution  system 
(including  LOGAIS)  in  terms  of  efficiency  and  effectiveness 
is  singularly  important.  Components  linked  together  as  a 
system  can  produce  an  end  result  which  is  greater  than  that 
possible  by  individual  performance.  Consequently,  all  com- 
ponents must  be  evaluated  on  the  basis  of  their  contribution 
to  the  system,  as  opposed  to  individual  performance. 

In  final  analysis  it  makes  very  little  difference 
whether  a firm  spends  more  or  less  dollars  for  an 
individual  component — warehousing  for  example — as 
long  as  the  overall  logistical  objectives  are  achieved 

as  the  lowest  total  cost  expenditure  [1;'17-'18].  j 

] 

Objective  j 

The  objective  of  this  research  effort  was  to  develop  ] 

a computer  model  for  minimiz,ing  LOGAIR  feeder  route  contract  ' 

costs.  A secondary  objective  was  to  vary  performance  para- 
meters in  order  to  analyze  trade-offs  between  cost  and  per-  I 

formance. 
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Chapter  2 


MODEL  DEVELOIMENT 


General 

IThe  mathematical  model  developed  for  LOGAIB  feeder 
routes  IS  a hiuary  linear  program.  The  model  was  designed 
to  find  a minimum  cost  route  structure  given  a fixed  number 
of  bases  and  various  performance  constraints.  The  model 
structure  allows  maximum  flexibility  for  perfoimiing  sensi- 
tivity analysis  over  a variety  of  constraint  requirements. 
The  characteristics  of  the  model  aire  explained  in  detail  as 
they  apply  to  the  model  derivative,  data  inputs,  and  vali- 
dation. The  software  package  used  to  solve  the  problem  is 
presented  in  SAUD  Report  EM-5627-PR,  and  is  also  available 
on  CREATE.^ 


-1 

A special  linear  programming  problem  whereby  all 
the  variables  are  restricted  to  two  values;  sero  or  one. 

A variable  is  used  to  indicate  whether  some  possible  action 
is  to  be  undertaken  (X='l ) or  not  (X=0)  (8:705). 

0 

“CREATE  is  an  acronym  for  the  GE /Honeywell  655 
computer  system  at  the  Air  Force  Institute  of  Technology. 
CREATE  stands  for  ^omputationsil  Resources  for  ^gineering 
^d  Simulation  ^Training  and  ^ucation. 
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Model  Derivation 


Previous  research  efforts  have  been  conducted  on  the 
general  vehicle  dispatching  problem  first  formulated  by 
Dantzig  and  Ramser  (2:309).  The  generail  problem  is  to 
satisfy  the  supply  demands  of  a set  of  customers  with  known 
locations  and  requirements  from  a depot  via  vehicles  of 
known  capacity.  The  objective  is  to  minimize  costs  of 
delivery  subject  to  constraints  of  customer  requirements, 
vehicle  capacities,  and  either  time  or  distance  restric- 
tions on  the  vehicle. 

A mildly  successful  evaluation  of  the  LOGAIR  freight 
delivery  system  within  the  context  of  the  vehicle  dispatching 
problem  was  initiated  as  early  as  1966.  Recent  efforts  to 
improve  the  LOGAIR  system  were  the  route  structure  model 
developed  by  Pailmatier  and  Prescott  (14),  and  the  subsequent 
improved  model  by  McPherson  and  O'Hara  (11).  Although  these 
models  considered  essentially  the  same  problem,  the  solution 
times  required  for  sensitivity  analyses  could  not  be  obtained 
from  the  branch  and  bound  algorithm  employed.  The  inef- 
ficiency of  the  branch  and  bound  technique,  when  used  with 
the  vehicle  scheduling  problem,  has  previously  been  noted  by 
Christofides  and  Eilon: 


^A  combined  heuristic  and  linear  programming  route 
generation  technique  was  proposed  by  the  RAND  Corporation 
in  October  1966  (5:4-11). 
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Clearly  the  computational  efficiency  of  the 
branch  and  hound  algorithm  when  applied  to  the 
vehicle  scheduling  problem  is  substantially 
reduced  when  compared  with  its  efficiency  in 
solving  an  equivalent  traveling  salesman  prob- 
lem [2:316]. 

The  inability  of  the  McPherson  and  O'Hara  model  to 
reach  an  optimal  solution  for  the  12  node  trunk  structure 
lends  credence  to  this  contention. 

Several  heuristic  models  were  investigated,  however, 
the  improbability  of  obtaining  consistent  near-optimal 
solutions  using  heuristics  eliminated  such  methodology  from 
extensive  consideration  (2:516-318).  Near  optimal  solutions 
would  be  necessary  to  conduct  the  desired  sensitivity  ana- 
lysis on  the  performance  parameters. 

A model  was  developed  by  this  research  team  which 
combined  the  fearures  of  optimality,  simplicity,  ease  of 
application,  amd  amenability  to  sensitivity  analysis. 

The  Model 

The  problem,  as  was  previously  stated,  is  to  mini- 
mize the  contract  costs  of  transporting  cargo  within  a 
LOGAIR  feeder  system.  The  present  feeder  systems  are  com- 
pletely airlift  dedicated  and  do  not  consider  other  methods 
of  shipment.  The  task  of  generating  a minimum  cost  feeder 
route  structure  involves  selectively  choosing  the  cheapest 
method  of  transporting  cargo  between  each  base  and  its 
parent  ALC,  providing  the  selected  mode  can  meet  the 
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requirements  of  time  and  capacity.  To  determine  the  totaJL 
minimum  cost,  consideration  must  he  given  to  all  possible 
routes  among  individual  bases  and  between  the  individual 
bases  and  the  ALC,  Also,  all  possible  modes  of  transporta- 
tion on  these  routes  must  be  considered.  Since  daily 
demands  of  cargo  at  each  base  are  known  and  the  cost  of 
each  mode  of  transportation  can  be  determined,  an  objec- 
tive function  was,  therefore,  formulated  to  minimize  the 
total  transportation  cost  of  a feeder  route. 

The  objective  function  is  subject  to  several  con- 
straints such  as  (1)  delivery  time,  (2)  minimum  vehicle 
visits,  (3)  out  and  back  requirements,  (4)  balance  con- 
straints, and  (5)  vehicle  capacities. 

Objective  function.  The  decision  variables  of  the  model  are 
defined  to  be  binary  variables  indicating  the  options  of 
traveling  from  one  location  (customer  base)  to  another  via 
a specified  vehicle.  A binary  variable  is  equal  to  one  if 
the  corresponding  route  and  vehicle  is  selected;  othei".\rise , 
the  variable  is  equal  to  zero.  The  sum  of  the  cost  coef- 
ficients of  the  non-zero  variables  in  the  obj'ective  func- 
tion then  determines  the  cost  of  that  feeder  route.  The 
objective  function  is  expressed  mathematically  as  follows: 
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Minimize : 


m 


2 2 1 

i=i  j=i  k=i 


i / d 


Where 


and 


C . = cost  to  travel  from,  base  i to  base 

^ G vehicle  k 


Z. 


i;jk  = 


, if  vehicle  k is  used  between 
base  i and  base  j 


^0,  otherwise 

i = departing  base  index 

^ = arriving  base  index 

k = number  of  the  vehicle  placed  in  the 
system  for  consideration 


Constraints.  The  first  set  of  constraints  are  the  delivery 
time  constraints.  These  are  included  to  allow  comparisons 
among  route  costs  for  various  maximum  delivery  time  criteria. 
Time  is  specified  in  terms  of  hours,  and  the  constraint  indi- 
cates the  maximum  time  allowable  to  deliver  goods  from  the 
hub  of  the  feeder  network  (the  ALC)  to  all  of  the  nodes 
(bases)  in  the  network.  That  is,  a vehicle  may  depart  an 
ALC  and  visit  one  or  more  bases,  but  it  must  return  to  the 
ALC  within  a maximum  specified  time.  Only  those  interbase 
routes  v;hose  decision  variables  are  not  zero  are  included  in 
the  summation.  The  general  delivery  time  constraint  is 
expressed  mathematically  as  follows: 
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Constraint  1 


n 


n n 

: 2 ^ ^ k = 1,2,3,... 

i=1  j=1 

Where  '^'iik  ^ time  in  hours  to  travel  from. 

^ base  i to  base  ;j  in  vehicle  k 

!T  = maximum  time  allowable  for  the  sum 
of  interbase  segments 

The  second  constraint  insures  that  each  base  will  be 
visited  oj  exactly  one  vehicle  since  the  sum  of  the  set  of 
decision  variables  for  all  vehicles  on  each  route  entering  a 
base  is  set  equal  to  one.  This  means  that  all  but  one  of  the 
decision  variables  must  equal  zero,  allowing  only  one  vehicle 
on  a particulajT  route  segment.  It  is  assumed  that  one 
vehicle  per  day,  either  aircraft  or  truck,  will  be  sufficient 
to  satisfy  the  daily  demand  of  any  base.  Constraint  two 
is  expressed  as  follows: 


n m 

2 ^ d=1,2,  3,  ...n 

i=1  k=1  i / d 


Constraint  2; 


The  out  and  back  constraint  insures  that  the  number 
of  vehicles  which  depairt  the  ALC  will  equal  the  number  which 
return.  The  constraint  is  expressed  mathematically  as  fol- 
lows: 


1^. 

This  constraint  is  partially  relaxed  at  the  ALC  to 
allow  all  solution  vehicles  to  return  to  the  ALC. 


1-4- 


n 

n 

n 

Constraint  3:  ^ 

I 

i 

i=1 

d=i 

i=1 

for  k = 1,  2,  3,  . . .m 

i 

"base  1 = ALC  H = large  nuxlier 

The  out  and  hack  constraint  for  a vehicle  "k"  insures 
that  if  that  vehicle  is  used  any’/nere  in  ahe  feeder  network 
it  will  he  forced  to  return  to  the  ALC.  In  conjunction 
with  Constraint  4,  this  will  also  insure  that  each  vehicle 
in  the  system  departs  the  ALC.  In  essence,  this  constraint 
forces  each  vehicle  out  of  the  ALC  and  then  hack  into  the 
ALC. 

The  fourth  set  of  constraints,  the  balance  con- 
straints, insure  that  a vehicle  which  enters  a base  also 
leaves  that  base.  These  constraints  are  written  as  follows: 
n n 

Constraint  4;  I ^ijk  I ^jlk  k = 1,2,  5,  ...n 

i=1  1=1  j =1,  2,  3,  . . eC 

i d 

If  vehicle  k travels  from  base  i to  base  j , it  must  also 
travel  from  base  j to  base  1.  In  other  words,  any  vehicles 
entering  base  j also  will  leave  base  j. 

The  capacity  constraints  are  similar  to  the  time 
constraints  in  that  they  can  provide  a comparison  of  route 
structures  controlled  by  vehicle  weight  and  volume 
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limitations.  These  constraints  will  insure  that  the  most 
restrictive  of  the  two  will  contribute  to  the  selection  of 
the  proper  decision  variables. 

n n 

Z Z for  k . -I,  2,  5,  ...  n. 

i=1  i / 

n n 

Z Z \ 

i=1  j=1 

V . , V.  = weight  and  volume  inputs  from  the 
^ ^ ALC  to  base  j 

V^,  V,  = weight  and  volume  capacities, 
respectively  of  vehicle  k 

The  weight  and  volume  constraints  provide  that  individual 
capacities  of  each  interbase  route-vehicle  combination  will 
be  summed  and  compared  to  the  maximum  allowable  vehicle 
capacity. 

Problem  size.  Problems  formulated  with  the  basic  model  will 
have  the  following  numbers  of  variables  and  constraints: 

Humber  of  vaniables  = I2]  x 2 x M (1) 

Humber  of  constraints  = 4-M  + H + M (H-1  ) (2) 

'^ere  H = total  number  of  bases  in  the  network,  and 
M = total  number  of  vehicles  in  the  network. 


Constraint  5: 


Where 
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ms 


■"IBBWi 


As  can  "be  seen  from  these  equations,  the  number  of 
variables  and  constraints  increases  exponentially  as  the 
numbers  of  bases  and  vehicles  increase. 

The  implicit  enumeration  algorithm  which  will  be 
used  to  solve  the  problem,  is  dimensioned  for  100  constraints 
and  150  variables.  Ko  further  increase  in  dimension  was  pos- 
sible because  of  the  memory  capacity  limits  of  the  computer 
system. 

Data  Collection 

The  LOGAIR  feeder  system  data  required  were  available 
at  ATTiC  Headquarters  and  also  within  current  Air  Force  publi- 
cations. The  data  were  categorized  as  follows:  (1)  customer 

demand,  (2)  vehicle  capacity,  (3)  vehicle  transit  time,  and 
(4)  direct  (vehicle)  operating  cost.  Each  of  these  cate- 
gories are  discussed  here  in  relation  to  information  sources, 
model  data  requirements,  and  data  transformations. 

Demand.  Daily  demands  for  every  major  Air  Force  base  are 
available  in  the  standard  data  base  of  the  AFLC  Control  Data 
Corporation  (CDC)  6400  CYBER  computer.^  The  base  demands 
represent  the  average  daily  tonnage  and  cubic  feet 


5 

'^The  demand  data  were  obtained  from  the  DoD  Material 
Distribution  Systems  Study  Group  (DODMDS),  courtesy  of 
AFLC/XRS  (16). 
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requirements  for  the  period  1 October  197^  through 
30  September  1975.  The  demand  data  require  no  transforma- 
tion, and  are  presented  in  Appendix  C.  Discussions  with 
personnel  at  the  AZLC  Directorate  of  Transportation,  LOGAIS 
and  Requirements  Branch  (LOTSL)  indicated  that  these  demand 
I figures  can  be  considered  accurate  estimators  of  demand  for 

I route  planning  purposes  (7;9). 

’{ 

;i  Capacity.  7ar*ious  vehicle  capacities  by  weight  and  volume 

are  required  for  establishing  the  capacity  constraints. 

Both  wei^t  and  volume  are  included  to  insure  that  the  most 
restrictive  feature  directly  influences  the  optimal  decision. 
Aircraft  capacities  were  obtained  from  the  LOGAIR  PT  77 
Plight  Schedules  and  Routing  Guide  (17:ii-iv).  Vaxious 
truch  capacities  were  obtained  from  the  Chief  Quality  Con- 
I trol  Inspector,  Transportation  Office,  at  Wright-Patterson 

APB.  Due  to  the  wide  variety  of  trucks  employed  nationwide, 
only  standard  capacity  vehicles  (i.e,,  the  most  common)  were 
considered.  Various  vehicle-  capacities  are  tabulated  in 
Appendix  C. 

Transit  time.  The  travel  time  for  any  vehicle  between  any 
two  feeder  nodes  is  equal  to  the  distance  between  the  nodes 
times  the  estimated  average  speed  of  the  vehicle  plus  the 
load  and  unload  time.  Air  miles  used  are  great  circle 
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straight  line  distances  (statute  miles)  obtained  from  cur- 
rent CAB  directives  (17).  Surface  miles  were  obtained  from 
AIM  177-1 35,  Official  Table  of  Distances  (20).  Distance 
matrices  for  the  air  and  surface  miles  and  the  average 
speeds  and  load  and  unload  times  for  each  vehicle  type  are 
shown  in  Appendix  C. 

Cost.  The  determination  of  costs  to  service  feeder  segments 
by  different  modes  of  transportation  is  different  for  air- 
craft and  surface  vehicles.  Aircraft  cost  data  were  ob- 
tained from  AlTiC,/LOTSL  and  include  the  basic  plane  cost  per 
mile,  fuel  subsidy,  adjustment  for  the  carrier's  net 
earnings  tax  and  landing  fees.  The  first  three  factors  were 
transformed  into  a total  per  mile  cost;  this  resultant  per 
mile  cost  times  air  distance  between  feeders  and  the  landing 
fee  produces  the  cost  per  route  segment  per  aircraft.  The 
cost  of  service  by  various  trucks  can  only  be  estimated  from 
existing  dedicated  truck  contracts  and  recent  contractual 
history,  as  a standard  rate  does  not  exist.  Each  service 
contract  is  negotiated  individually  by  the  Military  Traffic 
Management  Command  (MTMC)  (9).  Appendix  C lists  the  costs 
per  mile  for  both  aircraft  and  tracks. 

Model  validity.  Internal  validation  was  a continuing  process 
throughout  all  phases  of  the  model  development  and  testing. 
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Th.e  internal  validity  was  enhanced  by  continual  cross- 
checking of  mathematical  ma.nip'ulations  and  tabulations  for 
accuracy  and  correctness.  Assumptions  and  logic  flow  were 
continually  reviewed  for  appropriateness  particularly  in 
the  development  phase. 

Three  external  validation  phases  were  required  to 
insure  that  the  model  worked  properly  and  that  the  total 
software  package  was  convenient  to  use.  The  first  phase 
required  that  the  pilot  problem  be  manually  formulated  and 
put  into  the  EIP30C  program  to  insure  that  the  problem  could 
be  solved.  The  second  phase  required  the  validation  of  the 
matrix-generative  POSTRAJJ  program  which  was  written  to  facil- 
itate the  data  input  for  the  RIP50C  program.  The  third  and 
final  external  validation  phase  consisted  of  solving  for  an 
optimum  solution  for  the  present  LOGAIS  Route  and  perfor- 
ming a sensitivity  analysis.  This  phase  of  the  model  vali- 
dation is  covered  in  Chapter 

Sensitivity  analysis.  An  extensive  sensitivity  analysis  was 
conducted  to  complete  the  external  validation  and  to  demon- 
strate the  flexibility  of  the  model.  Minimum  cost  route 
structures  were  determined  for  a variety  of  input  conditions , 
generalized  as  follows; 

1.  Alternatives  to  the  present  policy  of  one-day 
delivery  to  each  customer  were  explored. 
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2.  Base  demands  were  vairied  to  include  priority  5 
big  load,  and  priority  5 big  and  small  load,  cargo. 

3.  The  type  of  vehicles  were  varied  in  an  attempt 
to  determine  a relationship  with  size  of  demands. 

Summary  of  Assumptions 

1 . The  data  obtained  from  APLC  Transportation 
Directorates  are  valid. 

2.  Demands  are  deterministic  (known  with  certainty). 
5.  Time  measures  are  linear  (11:23). 

4.  More  cargo  goes  into  a base  than  comes  out. 

5.  One  vehicle  per  day,  either  aircraft  or  truck, 
will  be  sufficient  to  satisfy  the  daily  demand  of  any  base. 


0 
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CtLapter  5 


PILOT  STULT 

General 

Although  the  basic  model  is  hi-directional  in  nature, 
the  problems  considered  here  were  formulated  to  permit  travel 
in  one  direction  only.  The  exception  to  this  rule  is  that 
the  arcs  connecting  the  bases  with  the  ALC  permit  travel  in 
both  directions.  This  exception  is  to  permit  a vehicle  to 
visit  any  base  in  the  system  at  either  the  start  or  the 
finish  of  a route,  depending  on  the  merits  of  the  individual 
arrangements . 

■The  one-way  travel  restriction  has  the  advantage  of 

'I 

reducing  the  nmnber  of  variables  in  the  model  and  thereby 
permits  a larger  problem  to  be  solved  with  the  given  com- 
puter limitations  (time  and  core  space).  If  each  arc  con- 
necting two  bases  is  allowed  to  have  traffic  flow  in  either 
2 

direction,  then  the  total  number  of  variables  in  the  prob- 
lem is  computed  by  use  of  Equation  1 , 


-1 

The  number  of  constraints  is  determined  by  the 
total  number  of  nodes  and  vehicles,  and  is  not  affected  by 
either  flow  selection. 

2 

It  IS  assumed  that  there  is  no  difference  in  cost 
between  traversing  an  individual  arc  in  one  direction  vice 
another.  The  shortest  connecting  route  between  several 
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Bj  contrast,  the  number  of  variables  in  the  problem 


foiTEulation,  given  the  one-way  restriction,  is: 


ITumber  of  variables 


(IT-1) 


X M 


(3) 


Por  comparison,  Table  1 and  2 summarize  the  number 
of  variables  required  for  various  combinations  of  nodes  and 
vehicles  for  both  one-way  and  two-way  travel.  The  area  below 
the  heavy  lines  indicates  the  area  ’^here  artifically  imposed 
computer  limitations  would  not  peimit  complete  problem  solu- 


tions. 


Pilot  Problem 

A logistics  network  which  consisted  of  three  cus- 
tomer bases,  one  ALC,  and  two  vehicles  was  fabricated  to 
assist  in  validating  the  model  and  developing  the  necessary 
software.  The  problem  was  constrained  so  as  to  permit  a 
solution  to  be  determined  either  by  inspection  or  by  an 
implicit  enumeration  algorithm.  Of  the  three  customer  bases 
(nodes  2,  3,  and  4),  two  of  them  (nodes  5 and  4)  were 
located  sufficiently  close  to  the  ALC  (node  1 ) to  permit 
access  by  surface  transportation.  The  other  base  (node  2) 
was  accessible  by  aircraft  only  because  of  the  time 


nodes  may  be  alfected  by  the  one-way  restrictions  however. 
The  optimality  of  the  solution  in  this  case  will  depend 
upon  the  numbering  system,  as  will  be  discussed  later. 


23 


Table  1 


Basic  Variables 
(Qae-Vay  Trarel  nodel) 


Humber  of 
Bases,  H 

Humber  of  Vehicles,  M 

2 

5 

4 

5 

6 

2 

4 

: 6 

8 

I 10 

: 12 

3 

10 

1 15 

I 

20 

I 

25 

30 

4 

18 

27 

56 

^5 

54 

5 

28 

42 

56 

70 

84 

6 

40 

60 

80 

100 

120 

7 

54 

81 

108 

1 

j 

155 

162 

8 

70 

105 

140  1 

175 

210 

9 

88 

152 

176 

220 

264 

10 

I 

108 

162 

216 

2?0 

324 

11 

130 

195 

260 

525 

590 

12 

154 

251 

508 

585 

462 

2^ 


Table  2 

Basic  Variables 
(Tv/o-Wa^’  Travel  Model) 


Nmnber  of 
Bases,  IT 

Number  of  Vebicles,  M 

1 

2 

3 

4 

r~ 

P 

6 

d. 

4 

6 

8 

10 

12 

3 

12 

18 

24 

30 

36 

4 

24 

36 

48 

60 

72 

5 

40 

60 

80 

ICO 

120 

6 

60 

90 

120 

150 

180 

7 

84 

126 

168 

210 

252 

8 

112 

168 

224 

280 

536 

Q 

</ 

144 

216 

288 

560 

442 

10 

180 

270 

360 

450 

540 

11 

220 

330 

440 

550  j 

660 

12 

264 

396 

528 

660 

792 
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constraint.  For  the  saJce  of  simplicity  only  two  vehicles 
were  considered  in  the  problem — an  aircraft  and  a truck. 

[The  values  selected  for  the  problem  variables  were  similar 
to  actual  data. 

Problem  Solution 

Jlanual  solution.  The  optimTim  solution  to  the  pilot  problem 
is  shown  in  Appendix  B.  It  can  be  readily  seen  that  there 
are  only  two  feasible  solutions  to  the  problem  because  of 
the  constraints.  The  optimum  solution  is  the  feasible 
solution  with  the  lowest  total  cost,  i.e.,  vehicle  1 (air- 
craft) going  from  base  1 to  base  2 and  returning  to  base  1, 
and  vehicle  2 going  from  base  1 to  base  3 to  base  h and 
returning  to  base  1.  The  total  cost  is  shown  to  be  $6,800. 

Comnuter  solutions.  In  order  to  solve  the  pilot  problem  with 
the  EIP30C  program,  it  was  necessary  to  set  up  the  problem 
in  the  following  standard  format; 

N 

Minimize ; I 

3=1 

N 

Subject  to:  b.  — 0 i = 1,2,3,.*.^I’ 

X k-j  — j <j 

j=1  X.  = 0 or  1 
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This  standard  format  required  that  all  the  con- 
straints of  the  basic  model  be  formulated  using  "greater 
than  or  equal  to"  signs  as  shown  in  Appendix  B.  The  same 
problem  is  also  presented  in  the  standard  format  required 
by  the  computer.  This  change  resulted  in  increasing  the 
number  of  constraints  over  that  which  was  required  by  the 
basic  model. 

Approximately  ei^t  hours  were  needed  to  manually 
set  up  the  constraints,  compute  the  vaniable  coefficients, 
and  type  the  matrices  into  an  accessible  file  for  the 
RIPpOC  program.  Additionally,  the  variables  had  to  be 
renumbered  consecutively,  i.e.,  ^®came  became 

Z2,  etc. 

The  pilot  problem  was  manually  put  into  the  EXP50C 
program  and  was  successfully  solved.  The  solution  con- 
firmed that  the  model  was  valid  and  appropriate  to  rhe 
problem  solution.  Different  solutions  were  obtained  as  the 
constraint  requirements  were  varied  indicating  the  appli- 
cability of  the  model  to  sensitivity  analysis. 

PORTRAIT  Program 

The  time  required  to  manually  set  up  the  pilot 
problem  was  considerable.  Likewise,  a considerable  amount 
of  time  was  expended  each  time  a vehicle  was  added  to  or 
deleted  from  the  system  during  the  sensitivity  analysis. 
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Unfortimately,  the  coefficients  of  the  variables  being 
eliminated  could  not  be  simply  dropped  from  the  matrix. 

The  coefficient  matrix  had  to  be  completely  redefined 
before  it  would  be  accepted  by  the  RIP30C  program. 

Although  this  effort  was  reasonable  for  the  validation 
process,  it  was  apparent  that  this  amount  of  set  up  time 
would  seriously  hajaper  any  practical  application  of  the 
model. 

In  order  to  alleviate  this  problem  two  FORTRAIT  pro- 
grams were  written.  One  program  was  written  to  develop 
distance  matrices  for  both  the  global  statute  miles  (air 
routes)  and  the  surface  statute  miles  between  the  various 
bases  and  between  each  base  and  the  ALC.  This  program 
accepts  the  distances  ia  free  format  and  generates  the  dis- 
tance matrices  in  the  format  required  by  the  second  FORTRAN' 
program. 

The  second  FORTRAN  program  was  written  to  generate 
the  input  data  file  for  use  with  the  RIP30C  program.  This 
program  automatically  computes  all  cost  coefficients  for 
the  objective  function  and  the  time  requirements  for  each 
vehicle  on  each  route  for  the  time  constraint  statements. 
This  information  is  then  stored  in  a matrix  of  the  proper 
format  for  use  by  the  RIP30C  program.  A complete  listing 
of  both  programs  is  presented  in  Appendix  A. 
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A complete  set  of  iastinictioiis  is  also  provided  in 


Appendix  A.  TMs  is  a step-by-step  outline  of  the  procediires 
which  must  be  followed  during  the  solution  process. 

The  solution  obtained  by  using  the  K)ETEAET  generated 
matrices  in  the  HIP30C  program  is  precisely  the  same  as  was 
obtained  by  using  the  EIP30C  program  alone  (manually  gener- 
ated and  inputed  data  matrices).  The  optimum  solution 
variables  are  printed  below  the  "Least  Z Before  Variable 
Change"  line  in  the  solution  output  as  shown  in  the  pilot 
problem  solution  in  Appendix  B,  The  solution  vaniables  ane 
listed  as  vairiable  indexes  and  refer  to  the  subscripts  shown 
in  the  data  input  for  the  pilot  problem.  For  example , 
variable  index  11  refers  to  variable  subscript  152,  which 
indicates  that  one  of  the  solution  variables  is  the  use  of 
vehicle  number  2 between  base  1 and  base  5.  The  time  and 
cost  of  each  individual  variable  is  also  shown. 

Summary.  The  FOETRAN  programs  reduced  the  time  for  genera- 
tion of  a solution  matrix  from  approximately  eight  hours 
to  about  fifteen  minutes.  The  most  beneficial  aspect  of 
these  FORTRAIT  programs,  however,  is  their  accuracy.  This 
accuracy  is  insured  through  the  simplicity  of  the  manner  in 
which  the  programs  receive  the  inputs  from  the  user. 

Figure  2 provides  a complete  overview  of  the  compu- 
tational programs  used  in  the  solution  process. 
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Model  Efficienc:; 


The  basic  model  utilizing  two-way  ancs  will  develop 
only  optimal  solutions  when  used  within  the  confines  of 
given  computer  limits.  This  is  because  all  possible  route 
combinations  are  considered  in  each  problem  formulation. 
Under  certain  conditions,  however,  the  baisic  model  will  also 
develop  infeasible  solutions  due  to  the  possibility  of  a 
vehicle  looping'^  between  two  customer  bases.  Problems  for- 
mulated with  the  one-way  model,  on  the  other  hand,  may  not 
produce  an  optimal  solution,  depending  upon  the  method  by 
which  bases  are  numbered  for  the  problem  formulation,  but 
will  always  develop  a feasible  solution. 

It  was  decided  to  pursue  the  objective  of  this 
research  effort  utilizing  the  basic  model  with  the  one-way 
restriction  for  two  reanons:  ("1)  the  greater  efficiency  of 

the  model  with  the  one-way  restriction,^  and  (2)  the  apparent 
power  of  the  model  to  consistently  produce  optimal  or  very 
near  optimal  solutions.  This  last  statement  is  supported  by 
the  results  obtained  from  the  example  problems. 


^The  solution  is  infeasible  because  it  does  not 
provide  a connecting  route  to  or  from  the  ALC. 


Computer  time  will  always  be  limited,  whether  by 
capacity  or  by  other  resources.  Since  the  objective  is  to 
develop  a mathematical  computer  model  to  minimize  the 
feeder  route  costs  it  is  felt  that  the  more  customer  bases 
and  vehicles  which  can  reasonably  be  considered  at  the  same 
time  will  ultimately  lead  to  the  lov;est  overall  costs. 
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In  most  cases  the  solutions  obtained  with  the  one-way 
restriction  appear*  to  be  optimal.  Sub-optimal  solutions 
which  were  obtained,  fell  into  two  general  types,  both  of 
which  were  attributed  to  the  method  of  numbering  the  bases 
in  the  network.  These  general  types  are  depicted  in  figures 
5 and  4-. 

The  sub-optimal  solution  depicted  in  Pigure  3 is 
the  least  consequential  of  the  two  types  because  it  is 
easily  spotted  by  inspection.  The  key  to  the  sub-optimality 
is  that  the  vehicle  crosses  its  own  path.  Whenever  this 
occurs,  the  bases  in  the  affected  route  should  be  re-numbered 
to  eliminate  the  crossed  path.  The  revised  solution  should 
then  be  optimal. 

One  of  the  solutions  depicted  in  Figure  f is  a 
sub-optimal  solution  of  the  second  type.  From  the  infor- 
mation presented,  it  is  impossible  to  determine  if  (a)  or 
(b)  is  the  better  solution  unless  the  distances  are  measured. 

Again  the  determination  of  a possible  sub-optimal  solution 
must  be  made  by  inspection.  If  it  is  determined  that  the 
possibility  exists,  then  the  bases  involved  should  be 
re-numbered  and  the  problem  re-solved.  Experience  with  this 
type  of  solution  indicates  that  in  order  to  insure  an  optimal  ] 

solution,  base  B,  which  is  common  to  both  triangles  DBG  and 
ABC,  should  be  re-numbered  at  least  three  different  times. 
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Base  B should  be  given  the  hipest  number  in  its  cluster, 
the  lowest  number  in  its  cluster,  and  then  the  number  it 
would  receive  if  the  bases  were  numbered  sequentially. 

Summary.  Although  the  problem  of  sub-optimality  does  not 
seriously  limit  the  use  of  the  model,  it  must  nonetheless 
be  considered.  There  are  two  possibilities  for  overcoming 
this  limitation:  (1)  incresise  computer  time  and  memory 

space  for  the  size  problem  necessary  to  solve  the  basic 
(two-way)  model,  or  (2)  use  a logical  numbering  system. 

The  first  possibility  should  not  be  completely  ruled  out. 
LOGAXR  contracts  are  awarded  on  a yearly  basis  and  as  such, 
may  be  of  sufficient  importance  to  warrant  a large  block  of 
computer  time  and  space  to  guanantee  optimality.  The  second 
method,  although  heuristic,  costs  very  little  to  implement. 

A logical  numbering  system  dictates  that  the  bases  be 
assigned  numbers  in  a continuous  cloclrArise  (or  counter- 
clockwise) manner.  If  a cluster  of  bases  does  not  form  a 
naturally  circular  path  or  other  natural  pattern,  the  prob- 
lem should  be  run  with  different  numbering  systems  and 
scrutinized  for  the  general  type  of  sub-optimal  solutions 
as  previously  mentioned. 

Summary  of  Limitations 

1,  The  model  requires  the  LOGAIE  feeder  netv;orks 
to  have  specified  depot  and  customer  bases , as  it  will  only 
be  optimal  for  a given  set  of  nodes. 


2.  Computer  core  space  available  for  the  solution. 


of  the  LOGAIR  model  was  limited  to  65K.  This  limitation 
subsequently  limited  the  size  of  problems  which  could  be 
solved  during  this  research  effort  to  a maocimum  of  150 
variables  and  100  constraints. 

3.  Infeasible  solutions  occurred  with  the  banic 
model  due  to  looping  between  banes.  This  problem  was 
eliminated  with  a modification  to  the  basic  model. 

4.  Data  inputs  into  the  iXDETRAIT  matrix  generators 
are  subject  to  the  following  limitations: 

A,  Format  statements  in  the  program  allow  for 
a maximum  cost  coefficient  of  9999.9999.  Any  figure  higher 
than  this  will  be  truncated  from  the  left  side.  Once  the 
matrix  is  generated,  however,  the  lost  digit(s)  can  be 
manually  replaced  prior  to  running  the  RIP30C  program. 

B.  During  the  compilation  of  the  cost  coef- 
ficients a distinction  is  made  between  air  and  surface  trans- 
portation. This  distinction  is  based  solely  upon  vehicle 
per  mile  costs.  Vehicles  with  a per  mile  cost  in  excess  of 
S2.00  have  an  additional  cost  of  3250.00  (landing  fee) 
added  to  the  product  of  leg  length  and  per  mile  costs. 

It  was  felt  that  $2.00  provided  a safe  distinction  between 
air  and  surface  transportation  costs  for  the  purpose  of 
including  the  landing  fee. 
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A caution  must  tlierefore  "be  o'bser'^-ed:  In  tlie  even- 


a surface  venicle's  per  mile  cost  exceeds  $2,00,  a landing 
fee  of  $250.00  must  be  manually  su'btracted  from  the  cost 
coefficients  prior  to  running  the  HIP50C  program,  or  the 
program  must  be  altered.  A similar,  hut  opposite,  problem 
will  occur  in  the  event  an  aircraft  per  mile  cost  is  less 
than  $2.00. 

C.  The  time  computed  for  a vehicle  on  a leg 
of  the  network  includes  an  additive  constant  to  allow  for 
loading  and  unloading  time,  refueling,  and  taxi  time  (air- 
craft). These  additive  constants,  which  are  "best"  esti- 
mates, are  one  hour  for  aircraft  and  three  hours  for  trucks. 
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Chapter  4 


SESDLTS  i 

< 

General 

LOGAIE  Houte  of  the  IT  76  LOGAIH  system  was 
chosen  for  analysis  to  complete  the  third  and  final  phase 
of  validation  of  the  model.  In  an  attempt  to  thoroughly 
test  the  versatility  of  the  model  the  results  of  model 
solutions  obtained  by  vazyiiiS  performance  level,  vehicle 
mix,  and  the  daily  demand  were  compared.  Experiments  with 
changes  which  would  expand  the  possibilities  of  the  model 
but  not  contribute  to  increased  demands  upon  the  computer 
system  were  conducted.  Such  manipulations  consisted  of 
(1)  doubling  per  mile  costs  and  capacity  of  vehicles 
entered  into  the  system,  and  (2)  changing  requirements  to 
indicate  vehicles  visiting  bases  every  other  day.  This 
first  manipilation  has  the  effect  of  doubling  the  number 
of  vehicles  being  considered  to  increase  capacity.  Because 
the  vehicles  are  considered  in  pairs  this  does  not  increase 
the  number  of  variables  in  the  problem  formulation.  This 
technique  woTild  be  of  assistance  in  formulating  a problem 

consisring  of  bases  with  very  large  demands  (ALC  network).  j 

The  second  manipulation  of  visiting  some  bases  every  other  j 
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da7  would  be  of  cLSsistance  ia  formulating  a problem  with 
bases  with  small  daily  demands. 

Phase  analyses.  The  number  of  computer  runs  that  would  be 
necessary  to  conduct  aJLl  of  the  sensitivity  analysis  of 
interest  is  enormous.  In  order  to  cover  the  widest  possible 
spectrum  in  the  time  available,  an  ordered  method  of  conduc- 
ting the  analysis  was  required.  The  plan  which  was  decided 
on  involved  a series  of  four  phases  of  analysis;  each  phase 
being  a building  block  for  the  rest. 

The  plan  resulted  in  a reduction  in  the  number  of 
options  which  needed  to  be  investigated.  It  was  obvious 
after  several  runs  that  problems  were  not  constrained  by 
time  limits  greater  than  15  hours  for  aircraft  or  84  hours 
for  surface  traffic.  Constraint  times  greater  than  these 
limits  were  subsequently  eliminated.  The  time  constraints 

for  surface  vehicles  used  in  the  majority  of  test  runs  were 

'1 

36,  60,  and  84  hours.  The  time  constraint  of  "^5  hours  for 
air  service  is  constant  for  all  problem  formulations  in 
this  chapter. 


'These  times  correspond  to  delivery  times  of  24,  48, 
and  72  hours  (one,  two,  and  three  days)  to  the  last  base  on 
any  routs  selected;  the  difference  of  12  hours  is  an  allo^v'- 
ance  for  the  time  required  to  travel  from  the  last  base  back 
to  the  AliC  with  three  hours  of  loading/unloading  time  at  rhe 
ALC.  This  last  leg  of  each  route,  though  it  is  not  con- 
sidered for  timing  purposes,  must  nonetheless  be  included 
for  its  overall  cost  contribution. 
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It  was  decided  to  limit  th.e  consideration  of  IP5 

cargo  in  the  demand  requirements  in  order  to  oetter  compare 

the  res'Jlts  of  the  computer  generated  solutions  with  the 

present  route  costs.  (The  LOGAIE.  system  presently  is 

designed  to  ship  TP1  and  2 cargo  on  a demand  oasis,  and 

2 

TP3  cargo  on  a space  aTailahle  oasis).  TP3  cargo  was 
considered  in  several  instances  to  determine  the  effect  of 
this  additional  demand  on  the  overall  costs. 

He  suits  of  .Analyses 

Phase  I initially  considered  aircraft  as  the  only 
vehicles,  and  then  considered  surface  vehicles  only.  The 
purpose  was  to  determine  which  vehicles  of  each  type  were 
the  most  cost  effective.  These  results  were  used  in  a later 
analysis  which  considered  the  air  and  surface  vehicles  in 
combination.  Consideration  was  directed  at  reducing  the 
number  and  type  of  vehicles  which  must  be  considered  for 
further  analysis. 


2 

The  entire  T?3  cargo  demand  was  considered,  as  wexl 
as  only  the  "big"  TP5  cargo.  "Big"  cargo  is  defined  as  that 
which  cannot  go  by  conventional  means  (UPS,  SEA,  Parcel  Pose, 
etc.)  (7). 

'"The  aircraft  considered  were  the  L100,  LI  88,  and 
DC-9.  Surface  vehicles  were  43,  ^0,  28,  and  18  LP  (Longi- 
tudinal Feet)  trucks.  % 
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Initially  all  types  oi  aircraft  and  trucks  were 
considered  in  th.e  formulations  of  tlie  problems.  Runs  were 
also  made  using  TP1  and  2 demands  only  and  TP1 , 2 , and  3 
demands.  Tbe  results  of  these  runs  indicate  that  certain 
vehicles  favored  certain  loads  due  to  compatibility  of 

4 

capacity  with  demand.  The  vehicles  whrch  were  in  the 
optimum  solutions  when  considering  only  TF1  and  2 cargo 
demands  were:  (1)  the  L188  aircraft,  (2)  the  4-0  LF  truck, 

(5)  the  28  LF  truck,  and  (4-)  the  18  LF  truck.  The  vehicles 
which  were  in  the  optimum  solutions  when  considering  TP1 , 

2,  and  3 cargo  demands  were;  (1)  the  L188  aircraft, 

(2)  the  DG-9  aircraft,  (3)  the  40  LF  truck,  and  (4)  the 
28  LF  truck.  Based  on  these  observations,  it  was  decided 
to  consider  only  these  vehicles  in  mosT;  of  the  succeeding 
anafysis. 

Phase  II  analysis  was  designed  to  investigate  the 
utility  of  the  different  options  available  in  the  RIP3CC 
program.  This  type  of  analysis  is  important  if  desired  to 
know  what  options  should  be  selected  under  what  prevailing 
circumstances  in  order  to  obtain  a solution  in  minimum  time. 


The  L0GA.IR  route  5^^  is  presently  sen/iced  by  the 
L188  aircraft  only.  Using  the  standard  costing  an  presented 
in  Appendix  C,  the  cost  of  this  operation,  as  determined  by 
this  program,  is  $3»‘500  per  day.  This  cost  figure  is  there- 
fore defined  as  the  reference  for  this  route. 


i 

I 

J 

I 

i 

I 

I 
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This  series  of  runs  investigated  the  use  of  the  program 
options,  both  singly  and  in  combination.  The  problem  was 
formulated  with  two,  three,  and  four  vehicles  in  order  to 
determine  the  relationships  between  the  number  of  problem 
variables  and  the  computer  run  time  for  each  option. 

Eesults  of  these  runs  are  included  in  Appendix  D.  The 
program  options  are  described  below: 

1.  L?  Stant  Option.  This  option  causes  the  EIP50C 
program  to  first  solve  for  a continuous  solution  by  linear 
programming  and  takes  the  initial  partial  solution  as  one 
determined  by  variables  which  have  a value  of  "zero  or  one." 
The  overall  effect  is  to  examine  all  roundings  of  the  con- 
tinuous solution  first. 

2.  Augmentation  Option.  This  option  specifies  a 
modification  of  Baias*  Rule  whereby  only  "free  variables 
corresponding  to  fractional  duai  variables  of  (L?s)  are 
considered  as  candidates  [6:5]." 

5.  Imbedded  LP  Option.  This  option  allows  for 
intemai  L?  solutions  at  selected  intervals  in  the  enumer- 

g 

ation  process  to  assist  and  accelerate  the  fathoming 
process  by  developing  surrogate  constraints. 


5 

^Egon  Baias'  algorithm  for  solving  binaiy  Imear 
programming  problems. 

solution  frequencies  of  one  and  eight  were 
recommended  by  the  user's  manual. 
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4.  Humber  of  Surrogate  Constraints.  Tbe  RIP50C 

user's  manual  recommends  a maximum  of  four  surrogate  con- 
9 

straojits. 

The  Phase  II  analysis  indicated  that  for  problems 
consisting  of  seven  bases  (such  as  the  LOGAIH  Route  5R) 
and  two  vehicles  (54  variables  and  46  constraints),  the 
imbedded  LP  (Linear  Programming)  option  (frequency  =8) 
and  the  augmentation  options  would  generally  provide  the 
fastest  solutions.  The  sane  seven  node  problem,  but  with 
three  vehicles  (81  variables  and  62  constraints)  yielded 
to  any  combination  of  options  which  included  the  LP  option. 
Once  the  problem  size  was  further  increased  to  four  vehicles 
(108  vsiriables  and  78  constraints),  feasible  solutions  could 
not  be  obtained  with  the  imbedded  LP  option  (either  fre- 
quency) taken  singly,  or  in  combination.  It  was  concluded 
that  the  imbedded  LP  option  and  the  augmentation  option, 
when  paired  together,  provided  the  fastest  solution  times 
for  problems  consisting  of  seven  nodes  and  three  or  less 
vehicles.  The  LP  start  option  appeared  to  provide  the  best 
solution  times  for  langer  problems.  Figures  5 and  6 show 
the  relationships  between  the  number  of  variables  and  run 


'^Descriotion  of  RIP3OC  ootions  from  user's  manual 

(6:2-3). 
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Figure  5 
-1 

Run  Time  Vs  Variables 
(LP  Start  Option  Only) 


LP  start  option  produced  fastes-1 
for  problems  greater  than  108  variables. 


solution  times 


44 


00 

Variables 


120 


Figure  6 

'I 

Iterations  Vs  Variables 
(LP  Start  and  ITo  Option  Only) 


'I 

Problems  with  greater  than  81  variables  could  be 
solved  vfith  LP  start  option  or  no  option  only.  Problems  with 
132  variables  required  approximately  35,000  iterations  (lar- 
gest problem  solved). 


times  and  iterations.  A complete  summary  of  these  findings 
is  presented  in  Appendix  E. 
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Phase  III  analysis  consisted  of  obtaining  Route 
solutions  with  vanious  mixes  of  air  and  surface  vehicles 
using  those  vehicles  indicated  by  the  Phase  I analysis. 

This  phase  also  considered  additional  vehicle  combinations , 
such  as  more  than  one  of  the  same  type  vehicle,  as  well  as 
other  combinations  which  did  not  surface  in  the  Phase  I 


,1 

I analysis. 


The  problems  which  were  formulated  and  solved  for 
this  phase  of  the  analysis  are  tabulated  in  Appendix  D. 

A comparison  of  costs  versus  performance  level  shown  in 
Pigures  7»  8,  and  9 is  particuleirly  significant  in  that  it 
demonstrates  the  accomplishment  of  the  secondary  objective 
of  this  research.  Data  points  for  these  curves  were  taken 
directly  from  the  corresponding  tabulations  in  Appendix  D. 
Considering  the  present  LOGAIP  Route  5^,  for  example,  the 
transportation  cost  for  a given  performance  level  can  be 
accurately  gauged.  A similar  assessment  of  inventory  costs 
versus  performance  level  (taken  in  concert  with  the  former) 
would  yield  an  accurate  estimate  of  a very  large  proportion 
of  the  total  distribution  costs  for  this  route. 

Phase  IV  consisted  of  computer  runs  which  were  made 
to  demonstrate  other  capabilities  of  the  model  such  as: 
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LOGAIR  Route 
TP1  & 2 Cargo  Orly 


10  20  30  ^0  50  60  70 

Maximum.  Support  Time 
(Hours) 


Pigure  7 

Daily  Cost  Vs  Performance  Level, 
TP  itkl. 


LOGAIR  Route  5R 

TP1 , 2 , & 5 (All ) Cargo 


10  20  30  40  50  60  70 

Maocimum  Support  Time 
(Hours ) 


Figure  9 

Daily  Cost  Vs  Performance  Level  , 
TP  1,2, §3 


1 . Doutlrag  the  daily  demands  as  if  deliveiy  would 
take  place  every  other  day, 

2.  Doubling  the  capacity  to  represent  doubling  the 
number  of  vehicles  on  a route,  and 

5.  Solviug  a large  scale  problem  (combining  existing 
feeder  routes). 

l!hese  computer  runs  were  made  primarily  to  experiment 
with  various  input  manipulations  and  also  further  demonstrate 
the  versatility  of  the  model.  Problem  formulations  and 
solutions  are  shown  in  Appendix  D.  The  computer  limitation, 
is  shovm  in  terms  of  variables  and  constraints  in  Appendix 
2.  The  area  below  the  heavy  line  indicates  problem  formula- 
tions not  solvable  due  to  the  limited  computer  storage. 

Several  computer  runs  were  also  made  in  an  attempt 
to  solve  a completely  different  netv/ork  problem.  This  prob- 
lem was  a combination  of  feeder  routes  and  5Q.  The  pur- 
pose was  to  demonstrate  the  capability  to  solve  larger  net- 
work problems,  i.e.,  to  optimize  a network  with  a signifi- 
cantly larger  number  of  bases  considered  simultaneously.  The 
air  and  ground  distance  matrices  and  the  daily  demands  of  the 
bases  in  route  5Q  are  listed  in  Appendix  C. 

Because  of  the  core  space  limit  the  combination 
problem  could  not  be  formulated  without  arbitrarily  desig- 
nating four  bases  as  transshipment  points  (Figure  10).  This 


proTslem,  with  three  vehicles  considered,  resulted  in  152 
variables  and  78  constraints.  The  formulation  of  this 
problem  and  the  computer  solution  is  listed  in  Appendix  D. 
The  cost  of  supplying  each  sub-feeder  base  from  its  trans- 
shipment point  was  added  to  the  minimum  objective  function 
value  at  optimum  routing.  The  resultant  cost  of  S9,3^5  is 
considerably  less  than  the  present  daily  cost  of  $11,160® 
to  operate  routes  5S  and  5Q  independently.  This  represents 
an  annual  savings  of  $662,475,  or  about  16  percent  of  the 
present  cost. 


Q 

This  figure  was  obtained  by  calculating  the  present 
cost  of  operating  route  5Q  by  use  of  the  cost  factors  shewn 
in  Appendix  C.  This  cost  was  added  to  the  present  cost  of 
operating  route  5R  (previously  calculated). 
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Chapter  5 


COUCLUSICH 


Siimnary 

The  model  developed  and  tested  lor  this  thesis 
provides  a method  of  minimizing  LOGAZS  feeder  route  trans- 
portation costs.  It  hais  proven  to  he  an  effective  and 
efficient  computerized  method  for  analyzing  trade-offs 
between  cost  and  various  performance  levels.  The  ease  v;ith 
which  the  model  maiy  be  used  allows  for  a general  comparison 
of  distribution  options  such  as  mixing  transportation  modes 
and  varying  performance  level  requirements.  For  example  if 
the  desired  performance  level  of  one  day  service  on  feeder 
route  is  extended  to  25.6  hours  and  tracks  are  used,  then 
the  daily  cost  of  this  route  can  be  reduced  from  $5 >600  to 
S5»S95  (Figures  11  and  12).  This  represents  a savings  of 
$622,325  pel*  year,  or  30  percent  of  the  present  cost  of 
operating  route  5S.  The  savings  possible  by  the  introduction 
of  surface  vehicles  and  the  relaxing  of  performance  con- 
straints do  not  imply  that  the  lower  cost  system  is  better 
than  the  present  system.  Inventory  costs  and  materials 
handling  costs  need  also  to  be  analyzed.  In  addition,  the 
impact  of  transportation  performance  on  mission  capability 


C3; 


Present  Route 


should  he  considered.  However,  an;7  additional  costs  nay  he 
insignificant  when  compared  to  the  estimated  savings. 

The  size  of  problems  which  were  solved  with  this 
model  are  limited  hy  the  computer  storage  capability. 

Problems  formulated  for  this  thesis  were  limited  to  150 

'I 

variables  and  100  constraints.  In  order  to  obtain  solutions 
for  problems  with  lairger  dimensions , it  was  necessary  to 

reduce  the  number  of  bases  by  the  use  of  transshipment  points  i 

! 

or  reduce  the  number  of  vehicles.  Such  methods  did  not 
guarantee  optimality  due  to  the  hueristics  involved  but  did 
provide  a technique  for  obtaining  hopefully  near-optimal  or 
optimal  solutions. 

! 

Recommendations 

Complete  optimization  of  the  LOGAIR  distribution 
system  appears  to  be  out  of  reach  because  of  the  many 
variables  involved.  However,  by  integration  of  individual 
optimal  feeder  routes  into  a total  aircraft  system,  sub- 
stantial. cost  reductions  may  be  realized. 

The  most  lucrative  area  for  further  study  of  the 
LOGAIR  distribution  system  appears  to  be  performance  require- 
ments. As  was  shown  in  this  study,  considerable  savings  can 


'I 

The  RIPJOC  program  must  be  re-dimensioned  to  take 
advantage  of  increased  computer  capability. 
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be  achieved  by  relaxing  the  performance  requirements  and 
allowing  expanded  usage  of  surface  transportation.  A 
feasibility  study  to  determine  what  effects  the  relaxing  of 
performance  requirements  will  have  on  combat  capability  and 
mission  effectiveness  is  desired. 

It  is  recommended  that  this  model  be  used  as  an  aid 
in  designing  future  feeder  routes  in  the  LOGAIB  system. 
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CALL  CSTC()SF(Li:iO,LA..\_\,CIJ<(I,J,-<)) 

2 "1  ? I 

2 '?  ) IF(LA.2C.l)  GO  TO  124 
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2)9  ! 125  COMTi;;UE 
21  !0  no  13.)  L=1 
21  M 3A(L)  =-l  . 
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P5CJ  PRIM  25,  .CIJ<(I,J,'<)  ‘ - 

251  3 210  CONTI  ROE 

2520  200  CONTINUE 

2530  DO  220  1=2, N 

2540  J=1 

2550  IJN=IJK+1 

2560  ;<X(I,  J,0=IJK 

2570  PRINT  25,T.J,<,IJ<,TIJ-:(I,J,'<),CIJ'<(I,J,<) 

2530  220  CONTINUE 
259  3 190  CONTINUE 
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3QQD  330  CONTINUE 
3010  IF(J.GE.2)  GO  TO  320 
3020  DO  350  1 = 2,  N’ 

3030  DO  360  :<=1,N 

3040  CALL  CSTRANT(LIN0,LC,  AC,AD,A5.;<C,  <X(I  .0) 

3950  360  CONTI  :.'UE 
3060  350  CONTI  NMJE 
3070  320  CONTINUE 
305OC  FOLLO/iS  NRITING'  OUT  ^ BACK  CONSTRAINTS  TO  SOL  VATR  I'/ 
3090  DO  370  :<=I.N 

3109  <C=’<C  + 1 

3110  DO  330  1=1 ,N 

3120  DO  390  J=1.N 
3130  IF(I.GE.J)  GO  TO  390  • 

3140  CALL  CSTRA.T(  LIN’D, LC,  AC  , AD , AE , '<C,  '<X  ( I .J.K) , -1  .0) 

3150  390  CONTINUE 
3160  330  CONTINUE 
3170  DO  430  1 = 2, N 
3180  J=I 

3 1 90  C A LL  CSTR  ANT ( LINO  . LC , AC , AD , AE . :<C , <X  ( I , J , K ) , 9 99 . C ) 

3299  400  CONTINUE 
3210  370  CONTINUE 

3220C  FOLLOWS  WRITL.’G  WEIGHT  CONSTRAINTS  TO  SOL  'LATRI.'': 

3230  no  410  <=1,,N 
324 T KC=UC+1 
3250  DO  42  0 1=1,  N 
3260  DO  439  J=  1 ,N 
3270  IF(I.JE.J)  00  TO  030 
3239  W=(-I  . )*VfJ(  J) 

3290  CALL  CST2 ANT( LINO , LC, AC. AD, AE.KC, <X( I , J,X)  , W) 

3300  430  CONTINUE 
331  0 42^''  CONTI.rUE 
3320  DO  440  1=2, N 
3330  J=1 

T 34 v/  = (“I  ' 

335  0 c A LL  C5T9  ANT ( L I NO . LC , AC , A r , AE , <C , NX ( I , J , 0 , 0 ) 

3360  440  CONTINUE 
3370  410  CONTINUE 

3330C  FOLLO'WS  WRITING  VOLUME  C(;NSTR  A I NTS  TO  SOL  ’OATH  IX 
339  0 DO  45 j -<=l,M 
34  00  •<C='<C+1 
3410  DO  46 O I=I,N 
3420  DO  470  J=1,N 
3430  IF( I .GE. J)  GO  TO  470 
344.0  V=(-l  . )-vVJ(  J) 

3450  CALL  CSTN  A' JT(  LINO,  LC,. AC.  AD.  AE,<C,  <X  ( I . J , , V ) 

34  60  4 70  CON’TINUE 
3470  469  CO.NTINUE 
3430  "'•()  45.}  I = 2,N' 

349}  J=I 


Pigu.re  14  (continued^ 


64 


* t / ti 


35'  1 

V=(-l .) 

*VJ{ J) 

351  I 

CALL  CS 

T2.a:<T(Li::o,lc,ac,  ad, r,  j,  o .'/) 

35.20 

4 30  go;; 

ciNOE 

3530 

450  CON 

TINU3 

3540C  FOLLOVJ3  /iOITINO  BALANCE  CONSTNAriTG  TO  SOL  " 

3550 

GO  ^90 

3550 

NO  5 0.1 

J=2,N 

3[37T  '<C=<C+I 
35  ’ 30  51  -5  i=i  ,r; 

359.,'  L=J 

35-^'’  IF(  I.GIZ.  J)  GO  TO  520 


351  0 

CALL  CST.'-'AN'TCLI 

NO,  LC, 

AC, 

« f'' 

■T.  . ♦ 

\E 

r/T> 

% \V./  f 

<X(  I , 

T 

sJ 

« 

X)  . 

1 .") 

3520 

520  IF(L.  GE.I) 

GO  TO 

51  0 

35  3.  J 

CALL  GSTNANTCLI 

;:o,  LC. 

A 

• 1 

AD, 

.\F 

VP 

XX(L, 

T 

0 , 

-1  .'*) 

354  ' 

51  ■ CONTINUE 

36  5 '-J 

CALL  cst-a\;;t(li 

NO, LC, 

AC, 

AD, 

.'.E 

,NC, 

XX  (L, 

1 

9 

'.) . 

365.) 

5-T1  CONTINUE 

3670 

490  CONTINIJF 

36=0 

DO  533!  N=l,.!‘ 

3590 

no  540  J=2,M 

37  J 

<c=;-'c+i 

371  ’ 

DO  551  I = r,N 

372  1 

L=J 

373’ 

I~(r.Gi^.J)  GO  T 

O'  56  V 

37-’.  1 

CALL  CST'-LANTCLI 

NO.LC. 

A 

i r* 
ni-  > 

A "t 

,’'C. 

XX (I  , 

J 

9 

<), 

-1  . . ) 

375  ■) 

55vi  INC  L.  OF.  I) 

GO  TO 

5" ./ 

375  1 

NALL  CSTNA.N'TCLI 

:'o,Lc, 

■'C, 

AD, 

AE 

,'<C, 

CXIL, 

T 

9 

V)  , 

1.0 

3 77-J 

55  ' CviNTMU- 

37  ’ J 

CALL  GST  ’ ANTd.I 

NO.LC, 

-VC , 

AD . 

AE 

,xc. 

'NX  ( L , 

1 

9 

1 . 0 

379  .' 

543  CONTI  ;U3 

3^'  ' 'J 

530  C..1NTINUF 

3‘=1  )C 

: F0LL0..3  F:;n’'AT 

STATE 

T*Q 

FO 

SO 

r ff  AT 

r 

1 

3':')  ; 

5 F(N^"AT(  14, 1 X, 

13, IV, 

13. 

13 

,1N. 

f 1 , r‘ 

« 

n.i 

X.I3. 

P.F.  3 1" 

ii,ix,i3,i;:,ii. 

1X,I1  , 

IN. 

1 1 , 

IX 

,6E'L 

OGA.n 

• 

I 

D 

= '4  .1 

15  FOrC'ATC  15,  IX 

, ’-(F  ?. 

4,1 

" ) ) 

305  ) 

2 5 F0  9 'ATC5N,  12 

, 1 2 , n 

,13 

N,I 

• 

IP', 

- 

CQ 

9 • 

.2,/) 

3-5,1 

35  F09'.:AT(  15, 1 X 

,4(2(  I 

3,  1 

X)F 

1 ) 

.3,1 

3 = 70 

^5  F09XAT(I5,1X 

,4(F5. 

1 , 1 

;<) ) 

55  FOO'.'T(V) 

3.-9  ' 

IF(  LC.~3,  1 ) GO 

TO  5 7'’ 

"JO  1 

NX)  5 = '*'  1 = 1,4 

391  3 

AF(I)=o  ). 

3 9 *10 

5 3"'  CONTIN'jF 

T 

J,/  J,  J 

..•■:iT,:(  1 1 ,35)Li:; 

C) . ( AC  ( Li  ) 

A P 

f ' 

( L 

E)  ,A 

’-KLE) 

« 

T 

E = 1 

.LN-I 

T 31  IV 
: ~ ♦ J > V 


T 1 ' * 

i « ' « 


2'~’-’r  1 or.:i-=u  G)+i  ’■ 

:-9  5 ' ."'ITFC  I I ,45)  LINO,  ( AF(  I ) , 1 = 1 .-i) 

29 6J  5 7'^  5T04 
39  7-'  NN’O 

393.'  5'J'^OOJTIN3  CSTC03F  ( LT  NM.  A . '.A  . 7.M 

3*7 ; or'ENFij.;  aa(  i ') 


Pigure  1A-  (continued; 


4'-'C''  iA(LA)=ZA 
4v;i  1 L,A=LA+! 

4rOA  TF(LA.LZ.4)  GO  TO  910 

4'-’3  ' WFITFC  n , 15)  Li;iO>,  ( AA(LZ)  , LZ  = 1 ,4) 

404  .1  L rK;=  LI.K)  + I ,) 

4053  15  FOOO-\T(  15,  1 X,4(F?.4,  1 X ) ) 

4'/:6v3  f_A=l 

407 J 910  lOFTUOH 

4^10  EHO 

4 ;90  S'nalUTriF  5C0FF(  LINO  , LB , ;^b  ,?n ) 

•^1  AO  nr'C';3i(M  aF(  I ■) ) 

4110  AB(L3)=Z:’ 

41 P A LB=L3  + 1 

4130  IF(L3.LF.4)  GO  TO  92'^ 

4140  ■.■i,PITF<  11,15)  Li;iO,(  AB(LZ)  , LZ=1  ,4) 

4 150  Lr:o=Li:.'o+!  ■•) 

4 150  15  FO  1MAT(I5,  iX,4(F9.4,  IX)) 

4170  l:"=i 
4130  920  OETUAM 
419.)  END 

4200  S'JZFOUTIilE  CSTOAjTC  LI  MO , LC , AC  , AF  , AE,  ZC,  ZD  , ZE) 
■t210  OI’^EKSIOM  AC(4)  , A,0(4)  ,AE(4) 

4 2^  ) AC(LC)=ZC 
423  3 A'MLC)=ZZ' 

424.)  AE(LC)=ZE 
4250  LC=LC+1 

4250  fF(LC.LE,4)  CO  TO  930 

4270  :';2ITi(  1 1 ,35)LrF0,  (AC(LZ) , AFC  LZ)  , AE( LZ)  ,LZ  = 1 ,4) 
423.'  Li’:o=Li;;o+i ’) 

4.290  35  F02 '.AT(I5,IX,4C2(I3,  lX)FI0.3,rO) 

43.. ‘0  LC=1 
431'  93'.'  .i’ET'J.FM 
43P  ! 'iF 


figure  04  (contijiued) 
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Tiiese  instructions  assume  th.at  the  tvro  FOKDEAIT  pro- 
grams, distance  and  solution  data  ma~rix  generators,  have 
been  input  to  the  users  iile. 

1.  DISTAKCS  J1AT5II:  GEIJEIlAirOR; 

SISTEn?  PORT  OLD  DIS^ 

(insure  that  dimension  statement  corresponds  to  correct 
# nodes) 

*RDR 

= Input  air  statute  and  surface  road  miles  for  the 
feeder  route  nodes  in  the  following  order; 


(Air  Miles)  (Surface  Miles) 

2.  Solution  Data  Matrix  Generator 
OTSTEM?  PORT  OLD  MATEST^ 

(Insure  that  dimension  statements  are  correct  and  tha 
Call  Attach  files  correspond  to  those  in  Distance  Mat 
generator  program) 

♦RUN 

= Input  data  in  the  order  requested.  When  all  inputs 
have  been  made,  the  program  will  print  the  variables, 
cross-reference  numbers  of  variables,  time  coefficients 


-1 

DIS  was  the  file  name  used  in  this  research  for  the 
PORTRAIT  distance  matrix  generator  program. 

2 

MATSST  was  the  file  name  used  for  the  data  matrix 
generator  for  test  problem  solution. 

Figure  Ip 

Instructions  for  Using  FORTRAN  Programs 
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and  cost  coefficients.  The  solution  data  matrix  v/ill  be  l 

saved  to  the  Gall  Attach  file  L0G5H.  i 

♦DOIilE 

SYSTEM?  PORT  OLD  LOG5S 
Check  and  correct  following: 

(A)  LINE  ECR  CORRECT  # VARIABLES,  CONSTRAINTS,  ' 

ETC . ' 

(B)  ?IIN  VISIT  "B"  COEEE  EOR  RETIIBNING  TO  ALO  EROM  ] 

1.0  TO  6.0  AND  COEEE  EOR  BASE  #5  EROM  1.0  TO 

5.0  (OR  MORE)  IE  TP5  DEMANDS  ARE  EICLDDED  IN  ' 

WT  & VOL  DEMANDS 

(C)  EXISTING  BOTTOM  LINE...CBECK  EOR  A COMPLETE 
GROUPING  OE  EOUR;  IE  NOT,  COMPLETE  LINS  'aETH 
0 0 0.  0 0 0.  DP  TO  EOUR 

(D)  LAST  LINE  SHOULD  BE  INSEPITED  AS  A GROUPING  OE 
EOUR  ZERO  GROUPS  AS  EOLLOWS;  5290  0 0 0.  0 0 0. 

0 0 0.  0 0 0. 

♦DONE 

5.  RIP500  Solutions 
SYSTEM?  CARD 

OLD  OR  NE7V  - OLD  AEIT . LIB/RIPPIJN , R 
READY 

*10$;IDENT: 

♦65S : LIMITS:? 5, 62K,^K  . . . (OR  UP  TO  LIMITS  NECESSAPM) 

♦1003:SSLECTA:  LOG5R 
♦RUN 


daily  weigh.!  aoid  volume  requirement 
in  tons  and  hundreds  of  cubic  feet, 


Distance  Unit  = Miles 


Figure  16 


for  each  node 
respectively 


lot  Problem 


f • 


i.  j k I.-.-  - - ^ V-'  * 1 

oYSTliM  OL  ' MA'IliJi' 

:^'i£A  ;y 

*K’U!i 

-iYr:^R  TO'rAL  3\ses  + ALC)  = i:rrEG.:R(N) 


ENTE3  TOTA'.  Yi’M  CLH5=  INTEuE3(  Y) 


E';TE3  4 VARIA3L.^S  I.i  INITIAL  PARTIAL  TOUT  I0i;=JPTI0.1  1 


li.TER  = J,IT  ;10  lYBEODtJ  LP;  1,1  P OPIEP.ilSE 


-iTTEP  = 'J.Ir'  COEFF  SIGNS  .lOPMAL;  1 , IF  SIGNS  PEyERSEG 


-NTcP  = J,FOP  INTERMEDIATE  OUTPUT  TO  APPEAR;  1,IF  OT '-P : i I S.i 


!;:.T^P  = 0,IF  a)ST  COEFF  ALL,  INTEG_IP;  I,  IF  OT'I-RNIS^ 


L'.T-P  = >%f)P  OTTER  AJG"ENTATION  PjLE 


L.T.:R  FPEjU-LIGY  OF  INTERMED I ATE  OUTPUT  = INTEGE 


',('7  J-^PER  SOUND  IF  ENONN 


L.T;R  = 0,0’?  'LAX  ^ SURROGATE  GlilSTR  A I NTS,  IF  XNCVNN 


EhT-R  = F^-TIJENCY  f^F  I’«FEODEJ  LP 


:NT.;P  ‘L\X  TIME  IN  SECONDS 
= 12' 

LiTuP  TME  NJMO ER  OF  T'ilS  ppO.iLE?'  RUN-  LOGAIR 

■ W 


5'igure  17  ( conLln'aed) 
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VA3IA3L-  r'JAX  > 

II  '3(;iJ7) 

Ji  J ; 

-.f'T'VrV  nr  r-r  V'.* 

.t. 

'.•4--VV 

I 2 I 

1 

*?  ’ ' 

-?2-j 

J < 

I J i 

2 

C. 

• 

DO 

] 1 

1 '■*  1 

\ : ? 

1 “t  1 

0 

i . ^ 

00 

2 3 i 

1 

*T 

3.  M 

32o 

.2  4 1 

•3 

3. 1 j 

• '* 

-K; 

3 4 1 

1 .23 

00 

2 1 1 

/ 

3. 

32 -J 

< ” 

3 1 1 

3 

1 .22 

■DO 

■'!  i ) 

^ 1 i 

') 

t .2  ! 

35 

i 2 3 

1 .3 

23.  -C 

1 ■ 

■ *1 

! 3 ? 

1 1 

5.dJ 

1 ’ 

1 

i 4 T 

1 2 

5 . 

1 

P 3 P 

i 3 

2 3.  2 

1 , ■ 

2 4 2 

14 

22.25 

1 .0 

\ \r  ■ 

3 2 

1 5 

•3.5,' 

1 '• 

2 1 2 

1;> 

23. 

i ; 

. 

3 1 2 

! 7 

0.0  ' 

1 ■' 

- 1 2 

1 3 

5 . 5i ' 

1 

J.  ' ‘ 
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li'O.) 

20  13 

0 

20  0 0 0.0 

120 

0 1 3 L0GMR99 

* i 

ICIO 

3230. 

00'0.j 

55(1.  O'KW 

350 

.1003 

325.1. 0303 

I ^2;i 

34  00. 

000  3 

530.0.1)  3250. 

0000 

550.  0000 

1 4j3(’ 

530. 

iMOn 

1 

1 00 

. 300.) 

1 00.  3030 

1 »4,) 
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1 050. 1 Wl 

|O0 

.J0O.) 
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1 00. 
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1 J6.; 

12. 

) 

24 . 0000 
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-1 .0000 
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-1. 
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-1.0000 

6 

, 0030 
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1 . .1.100 

3 

• 

a. 
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1 7. 
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21 .0000 

30 

. 0003 

23. 0000 
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2 
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3 
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3 
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6 

14 
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6 

15 
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7 

7 

-1.000 

7 
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7 
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Pig-ire  17  C continued) 
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Table  ^ 

Pilot  Problem  Solution 


Possible 

Solutions 


Aircraft 


Truck 


1-2- 3-4-1 

7600 

1 -2-3-1 

7050 

♦ 

1-2-1 

*♦6500 

* 

1 -2-4-1 

7200 

1-4-1 

O 

o 

200 

1 -3-4-1 

1650 

**300 

1-3-1 

1100 

o 

o 

*Hot  feasible  because  of  time 

♦♦Optimal  combination  of  aircraft  and  truck 
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Table  5 


Pilot  Problem  Pormulation 
(Basic  Model) 


Let:  X.  = 1 , if  from  station  i to  station  j via  vehicle 


k 0 , if  not 


t-,-v  = ■oioelino  time,  including  onload/offload,  from 

* * - - 


station  i to  station  j via  vehicle  k 


w . = dailj  weight  demand  at  station  j 
V . = daily  volume  demand  at  station  g 

<J 

T,  = mazcimum  pipeline  time  via  vehicle  k 
= maximum  weight  capacity  vehicle  k 
7,  = maximum  vol\ime  capacity  vehicle  k 

Objective  fimction: 

Minimize:  3250X^21  + 31  + 

3400^2^^  + 55CZ^^  + 3250X2^^  + 550X_^^  + 

+ i000X^Q2  + 

1000X2^2  ^ 1 OOX^/,  ^ + 1C00Xq^p+ 

1 OOX 2 ^ 2 

Subject  to: 

Time  Constraints 

3.00X^21  + 1.20X^,^  + 1.20X^_^,,  + 3.00X23^  + 

2«10X2^x|  + '1.20X;j^^  + 3«00Xo.i^  + ^.20X-^^  + 
1.20X^^,  ^ 12 


?8 


Table  5 (continued) 


28.00X, 


22  5»5CX^2^2  + 


28COZ252  + 29.25X242  5. 50x^2  + 

28.00X2^  p + 5.5CX-^2  + 5.50X_^^^  — 24 


MUT  VISIT  COHSTRAHTTS 
121 


X/^  P^  + Xxi 


= 1 


^151  ^231  ^132 


^2^2  “ 


X, 


+ Xo;,xi  + X^,,^  + X^  ,,  O + X^ 

X-,,.^  = 1 


141  ^ ^241 
342 


^541  ^142  ■*■  ^242  ^ 


X21  ^ + Xj^  ^ ^1  + X2x|  2 + ^ 2 '*’ 

X412  ^1^6 


OUT  AITI)  BACK  COUSTEAHTTS 


X121  + x^,^ 

+ ^-141 

+ 223^ 

^241  '*'  'I 

+ ^3^^ 

+ 

999  (X2^.  + 

^311 

^41l) 

X122  + 2^52 

+ ^142 

+ 2^232 

^242  '*’  ^21 2 

+ 21j^2 

+ 2:^'!  2 

999  (X2^2 

^31 2 '*' 

2^412^ 

COUSTRAHTTS 

^121  ^231 

+ ^241 

+ X2^^ 

^1  31  ^251 

" ^31^ 

+ 2;^^ 

X^4^  + ^94-1 

+ 2;^^ 

= x^.. 

+ X 


341 


y-rc. 
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i 

i 


CAPACITY  CONSTRAHTTS 

^131  ^^141 

+ ^^241  ~ 


+ 2X 


251 


+ 


^0X^22  ■*■  ^^^52  "*" 


^^^121  ^ ^151  ^^141  '*’  ^251  ^ 

■*■  ^^241  ~ 

^^■^122  "*"  "^152  ^^142  "*"  ^232  '*’ 

5^^2  '*’  ^^242  “ 
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Table  6 


I 


! 


I 


Pilot  Problem  Poraulation 
( Traasi o rmed ) 

MirLimize : 

5250X^21  + 550X^3^  + 550X.^^  + ?25CX23^  + 
,^00X2^^  + + 5'^50^2'^-^  5^2X3^^  + 

550X^.,.,  + 1C00X^2o  + ^ 42 

1000X^^2  1050X^^2  + '1C0X3^P  + 1000X2^5  + 

^ 00X3^  2 o 


Subject  to: 


-3.0CX^2'1  - ^-20X^3^  - 1.20X^^^  - 3.00X^^^ 
-5.10X2^.  - 1.20X3^^  - 3.00X2.^^  - 1.20X3^3 
-1.20X  + 12  2:0 


-28.00X^p2  ~ 5.50X^3P  - 5.50X^^2 

-23.00X232  " 29.25X2^2  “ 5.50X3^2 
-28.00X2^2  “ 5.50X3,^2  - 5.50X^^2  + 2^  i 0 


^121  ■*■  ^122  “ ^ ^ 

-X^  21  ”^122'*’  ''  ” 

^131  ^231  ^132  ^232  - - 0 

-^'1^-1  - ^231  ""  ^132 


+ 1 i 0 


^'I4'l  ■*■  ^41  '*’  ^341  ''”  ^142  '*'  ^242  '*’ 


“TM-P 


0 
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Table  6 ( continued) 

“ ^341  “ ^142  " ^242  “ 

^542  ^ ^ 


■*■^211  ■*■  ^411  ■*■  ^212  ■*■  ^312 

"*'^4'12  “ ^ 

"^11  “ ^311  " ^411  " ^2^12  ~ ^512 


1 


f 

I 

i 

Table  6 (continued)  ; 

^32  ^252  “ ^512  " ^542  “ ° 

“^132  “ ^32  ^312  ^3^2  ~ ° 

^142  ■*■  ^242  ■*■  ^342  “ ^411  ~ ° 


“^142  " 

2242  2; 

542  ^411 

^ 0 

-10X^21 

- 2^131 

^^141 

2X231 

“ ^^241 

17  ^ 0 

-102^22 

- 2^132 

- ^^,42  “ 

2X252 

-^^342  ■ 

- ^^242 

21  i 0 

-152,21 

- ^,3, 

“ 52,^,  - 

2X231 

-5X54,  ■ 

- 3X24,  + 

30  ^ 0 

-132,22 

- 22,32 

- 52,^2  - 

2X232 

CM 

I 

“52242  + 

23  ~ 0 

S3 


I 

5 


i. 

1 


Table  7 

Pilot  Problem  Solution 


LOGTEST 

IMPLICIT  EErmnaA-TiOEr  complete  total  time=5173.59^ 

LEAST  Z BEPOBE  VABIABLE  CHANGE  = 6.80000000E  05 
1000007000  000  15 

0 0 18 

NO.  PEASIBLE  SOLUTIONS  ^ 

ZS  GE  ZBAR  0 TIMES 
CONSTEAINT  INPEASIBLE  5 TIMES 

AUOffiNTATION  IMPOSSIBLE  0 TIMES 
AUGMENTATION  POSSIBLE  8 TIMES 

INTEGER  DUALS  0 TIMES 
NO.  OP  ROUNDED  INT.  DUALS  0 

LP  PATHOMED  0 TIMES 
LP  CALLED  0 TIMES 
NO.  ITERATIONS  17 

LAST  PEASIBLE  SOLUTION  AT8S1 5.875  SECONDS 
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APPENDIX  C 

LOGAIR  PROBLEM  SUPPORT  DATA 
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^Average  daily  demands  less  than  .005  tons  or  .5  cubic  feet  are  indicated  by 

(SOURCE  OF  DATA:  DoD  HATEEIAL  DISTRIBUTION  SISTEM  STUDY  GROUP  (DODMDS),  Courtesy  of 
AFLC/XRS  (16). 


Feeder  Route  5Q 


Table  9 


Daily  Requirements; 


1 


Feeder  Bases^  Priority  1 & 2 

Tons  Cu  Ft 


Tyndall  AFB  FL  (PAM) 

1.16 

200 

MacDill  AFB  FL  (MCF) 

1.06 

175 

Patrick  AFB  FL  (GOF) 

.21 

57 

Homestead  AFB  FL  (HST) 

1.45 

310 

Key  West  HAS  FL  (HQK) 

.35 

66 

Jacksonville  HAS  FL  (HIP) 

.46 

128 

Moody  AFB  GA  (VAD) 

.14 

20 

'I 

Priority  5 cargo  not  considered  for  Feeder  Route  5Q 

2 

All  Offline  Support  Bases  not  considered  significant. 
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Tal3le  12 

Truck  Transportation  Costs 


Truck 

Estimated  Contrac 

p 

Type 

Per  Mile  Cost 

45IP  Trailer 

$1.50 

40LP  Trailer^ 

1.25 

28LP  Trail er/Yan 

1.05 

18LP  Van 

.95 

Most  common  carrier  vehicle  utilized  (9). 

2 

Estimated  contract  per  mile  costs  ane  based  upon  a 
review  of  existing  USAE  dedicated  contract  truck  routes — all 
routes  negotiated  by  MTMC.  Round  trip  contracts  which  were 
reviewed  included;  (1)  Tinker  APB  - Altus  APB,  (2)  McClellan 
APB  - Castle  APB,  (3)  NAS  Dallas  - Kelly  APB,  Tinker  APB  - 
Kelly  APB  - Dyess  APB,  (5)  Wright -Patters on  APB  - Rickenbacker 
APB,  (6)  Corpus  Christi  - Kelly  APB,  and  (7)  NAS  Dallas  - 
Kelly  APB. 


92 


g 

o 

CA 

LA 

CN 

LA 

00 

A 

CD 

00 

CD 

A 

LA 

3 

o 

LA 

CA 

4" 

CM 

4" 

A 

A 

O 

CA 

CO  O 

> 

o- 

CD 

CA 

CM 

CM 

CM 

4- 

A 

CD 

CN 

LA 

A 

P=H 

CD 

c- 

CM 

00 

O 

<r~ 

00 

r" 

CD 

(A 

A- 

A 

o 

O 

CD 

00 

C“ 

[A- 

-u 

4" 

CO 

CA 

CD 

CD 

O 1 

o 

CA 

00 

cD 

4- 

4 

C\J 

A 

00 

(A 

A 

O 

CO 

a^ 

LA 

00 

CD 

4" 

A- 

CD 

4- 

CM 

CD 

CO 

CM 

A- 

r“ 

v 

04 

O 

O 

u^ 

O 

1 1 

[>- 

CM 

CA 

[>- 

[>. 

■4 

V” 

LA 

A 

OJ 

C“ 

V 

g 

CD 

lA 

I>- 

CD 

A' 

K^ 

00 

CN 

3 

CD 

rA 

00 

c- 

2) 

O 

li^ 

T~ 

C 

o 

1 

1 1 

PQ 

CD 

CM 

rA 

LA 

LA 

CA 

CO 

■>* 

*T“ 

P=! 

CD 

(A 

LA 

rA 

A 

o 

M 

[>. 

C\J 

LA 

00 

00 

LA 

o 

1 

1 

I 1 

m 

LT\ 

V" 

lA 

00 

O 

A- 

A 

Ph 

CM 

O 

CM 

O 

CA 

CD 

H 

LA 

A 

4- 

00 

A 

r* 

O 

00 

O 

1 

1 

1 

1 1 

k 

CCJ 

00 

[>. 

4" 

A 

lA 

A 

LA 

EH 

CA 

CO 

4" 

CD 

[> 

LA 

£3 

a^ 

CD 

O 

LA 

CA 

CA 

o 

1 

1 

1 

1 

I 1 

w 

LA 

O 

CO 

CD 

LA 

A 

c- 

CD 

CD 

rA 

LA 

CD 

o 

CTc 

LA 

O 

O 

I 

1 

1 

1 

1 

1 1 

^A 

00 

CO 

LA 

4" 

CM 

00 

00 

CM 

3 

O 

hA 

(Jc 

CO 

O 

1 

1 

1 

1 

1 

1 

1 1 

KN 

CD 

*x 

CM 

CO 

O 

o 

CD 

CM 

1>- 

CA 

CD 

o 

1 

1 

t 

1 

1 

1 

1 

1 1 

> 

rA 

LA 

-1. 

-s 

A 

CM 

o 

H 

LA 

c>- 

O 

1 

1 

1 

1 

1 

1 

1 

1 

1 1 

pq 

•d" 

CM 

s 

t>- 

roi 

a 

o 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 1 

3 

CD 

CM 

& 

•ch 

fc 

O 

1 

1 

1 

1 

1 

I 

1 

1 

1 

1 

1 1 

\D 

H CO  a EH  CM 
H 0(  O CO  H 


Air  Statute  Miles 


^rCXDU^^OU^^':-OC^JC\JL^^^ 

vLrsLfNtAoji:-c\jij-^Lr\vD-:l'C^ 


C^Jr-C04•a^vDOOC^JN^OJv- 
LTNOj^-ajooajojcT'LrNvoroivDo  I 
N^C0[>-LrNrcsroiv-^'r-C^C0C\j 


[^-o^oa^^c^'r^<^ovDa^ 
vO^-'T-C^O^rOvOONCMOO  I 1 
L^^cy^oOLf^Lr^4-c^J^K^c^>a^ 


o-oD  ^ 

O0000v-^0^<^^<^<MC^JO!  I I 
LfN^  V L^^c^c3^c^^<^ 


cr 

LA 

[>- 

LD 

<T. 

<T> 

00 

cu 

4- 

LA 

0) 

CM 

1>- 

O 

CM  O 1 1 1 1 

pq 

N^ 

^<^ 

4* 

4* 

[>- 

C3N 

LD 

PL| 

V“ 

o 

r* 

hA 

O 

A1 

LO 

o 

ON 

LA 

S^ 

fA 

O- 

LA 

O I 1 1 1 1 

OvJ 

X) 

LO 

C\J 

CM 

LA 

5^ 

CM 

00 

CM  4: 

3 

V 

O 

C\J 

O X) 

PM 

CM 

lA 

CM 

CQ 

CO 

o- 

CM 

LA 

CM 

-s 

LD 

4-  O 

> 

CM 

3- 

LA 

r- 

T-  O I I I 
lA 


O I I I I I 


^ g ^ 

pq  (t 


C3N  IX)  OJ  C\J  CO  O 

C\JlA4-CMO^v-OI  I I I I 1 
LA  0\  CX)  X)  LA  ^ OJ 


lD  CO  CA  LA  rA  o 
NA4-^AJ04-0I  I I I 1 I I 
KN  iX)  4-  fA  OJ 


O ^A  4- 

lALDC\JOIIIIIIIlll 
^A  OJ 


LA  CA 

S^  4-  O 1 I I I 1 I I 

LA  Al 


APPENDIX  D 
COMPUTER  SOLUTIONS 


95 


Feeder  Route  5R  Analysis — TP  1,  2,  and  5 Big  Cargo 


0)  o 
E^ 

>■ 

+3 

u 

o 

ft  <U 
ft  iH 

0) 

> 


03 

03 

••  -P 

P!  a 
o ri  o 

0+3  0) 
H ftfH 
O O O 

oa  -H 

0) 

> 


iAr~  ^ lA, 

J- J)  J-  ki)  !>•  r*  ^ ^ 

llll  llllll  I I 

U?  OJ  4-  CM  OJ  in  CN^t  ^ CvJ  CM  CD  CM  U)  CM  CO 

III  IIIIIIIIIIIIIII 

PWftWWWWOftPftO<3)W<JQ«s!0 


rcifov^  V I 
111' 
CM  CM  CD  CN; 

llll 


Jr- Jv- 

I 1 

CM  CD  CM  CD 

llll 


98 


Combined  Feeder  Routes  5R  and  5Q  TP  1 and  2 Cargo 


RIP30C  Program  Limits  (Basic  Variables 
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Options  Analysis  Summaiy  for  5^  Variables  and  46  Constraints 
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